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In a recent Letter jlj, an attempt is made to justify 
and to provide a further support to earlier erroneous 
claims, see e. g. in Q, that particle distributions 
or, equivalently, probability densities can become non- 
normalizable in the case of anomalous Levy walk dif¬ 
fusion, and this anomaly can reflect a physical reality. 
In this Comment, I show that such claims cannot be 
justified and supported. 

To begin with, any probability density function or 
PDF must be non-negative and normalized. Both a 
common sense and the postulates of probability the¬ 
ory fix these important properties, once and forever, 
see e.g. Ch. 2 in 0]. In Ref. [lj], Eq. (8) is derived 
from the correct result in Eq. (2) therein by taking 
an asymptotic limit t — > oo. The PDF in Eq. (2) is 
obviously normalized for any particular model within 
the Levy walk model considered. The authors seem 
perfectly understand this, see their Eq. (3). However, 
already this fact implies that Eq. (8) cannot be cor¬ 
rect being applied to the whole range of x variations. 
Indeed, Eq. (8) describes a part of the distribution, 
for a sufficiently large x, however, not the entire dis¬ 
tribution. In this respect, Eq. (1) in [I| also describes 
only a tail of the time distribution ip{r). Clearly, the 
distribution in Eq. (1) in IdJ is non-normalizable, if 
by mistake to extend it to small r —> 0. However, to 
claim this for possible were absurd. An example of 
PDF can be readily given. It is Pareto distribution, 
ip(r) = a/[T 0 (l+T/To) 1+ “], where a constant r 0 sepa¬ 
rates small, t To, and large, r > ro, time intervals 
between the scattering events. By the same token, the 
claim that the distribution PA(x,t) in their Eq. (8) 
is not normalizable because PA(x,t ) ~ |a;| _1_a in the 
vicinity of x —> 0 is obviously wrong. The correct dis¬ 
tribution does not behave in this way for |x| < vqTq. 

After deriving incorrect result in Eq. (8), an at¬ 
tempt to further substantiate it is made in [l[ by using 
the concept of infinite covariant density (ICD). Below 
I explain why this attempt fails. For this, let us clarify 
the origin of ICD in their Eq. (9). It is not a PDF, 
but a function concocted to be non-normalizable for 
a > 1. Indeed, let us change random variable from x 
to v a = x/t a , where t is just a parameter, and not a 
random variable. Then, the corresponding PDF is 

Pj; a) (v a ,t) = J 5{v a - x/t a )P(x,t)dx (1) 
= t a P(x = v a t a ,t) . 

Most obviously, it is normalized for any t, 


fPf; a \v a ,t)dv a = 1. Hence, P^\v a ) := 
Hindoo Pv° \v a , t) is also normalized. At the first 
look, our P^\v a ) may seem to be their ICD. How¬ 
ever, this is not so. The function I c d(n) in Eq. 
(9) in [l[ is not a PDF of v = V\ , which would 
be linp^oo (ui, t) = lim^oo tP{x = vt,t) = 
Poo\vi), but a special construct. It can be consid¬ 
ered as the t —> oo limit of the function I c d(v,t) := 
t a ~ 1 Pv 1 \vi,t). Its normalization, 


J I cd (v,t)dv = t a \ 


( 2 ) 


indeed diverges for a > 1 in the limit t —> oo. Of 
course, at any fixed and finite observation time t, the 
(generalized) moment averages 


<M 9 > 


J \v\ q I C d(y,t)dv 
f I c d(v, t)dv 

J \v\ q Pf; 1) (v,t)dv 


= t 


1—a 


\v\ q I cd (v,t)dv 


( 3 ) 


can be found either with I cd (v,t), or with PDF 
Pv l \v,t). However, it is incorrect to use for this pur¬ 
pose I c d(n) obtained in the formal limit of infinite t, 
as the authors do in their Eq. (13). It is generally, e.g. 
for q < a, wrong. The limit t —>• oo must be taken in 
our Eq. ©. Importantly, stochastic numerics can be 
done only at finite t and ro. In this respect, which con¬ 
crete PDF was used in the numerics done in [l] 
is dim. Most obviously, any decent experiment, either 
real or numerical, done at finite t will yield I c d(v,t), 
which is perfectly normalizable, and not I cd {v )• Hence 
ICD cannot correspond to a physical reality, contrary 
to the claims in [lj]. 
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